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We developed a density matrix renormalization-group technique to study quantum Hall fractions 
of fast rotating bosons. In this paper we present a discussion of the method together with the results 
which we obtain in three distinct cases of the narrow channel, cylinder and spherical geometries. In 
the narrow channel case, which is relevant to anisotropic confining traps in the limit of extremely 
fast rotation, we find a series of zero-temperature phase transitions in the strongly interacting regime 
as a function of the interaction strength between bosons. We compute energies and density profiles 
for different filling fractions on a cylinder and compare the convergence rates of the method in the 
cylinder and a sphere geometries. 
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I. INTRODUCTION 

Density matrix renormalization-group (DMRG) meth- 
ods play an important role in numerical studies of one- 
dimensional quantum lattice models allowing to access 
their low-energy properties with a very high precisiouji 
Recently, these powerful techniques have also been suc- 
cessfully applied to nonequilibrium problems in strongly- 
interacting quantum systems such as dynamics of quan- 
tum quenches^ and tunneling of electrons through quan- 
tum point contacts.'^ A considerable theoretical effort 
has been put into development of similar methods for 
two-dimensional quantum lattice models, see for exam- 
ple Refs;-^ and** These new numerical techniques make 
use of the recent advances in understanding entangle- 
ment properties of many-particle quantum systems in 
finding the optimal basis for the representation of the 
ground state and low-lying excitations. This is usually 
implemented by using matrix-product states or tensor 
networksi^ DMRG methods have also been used in quan- 
tum chemistry and nuclear physics 

An interesting application of DMRG was developed by 
Shibata and Yoshioka^ and more recently by Feiguin et 
al. in the studies of interacting electrons in the frac- 
tional quantum Hall effect (FQHE) regime. In the paper 
in Refji, the authors calculated ground state energies, 
gaps and correlation functions of electrons at v = 1/3 and 
I' = 5/2 filling fractions in the spherical geometry for sys- 
tems with larger number of particles than it was possible 
to study with exact diagonalization. These calculations 
showed a potential of the method in finite-size numerical 
studies of quantum Hall systems. A standard approach in 
FQHE is exact diagonalization which has been very fruit- 
ful in identifying states of interacting electrons with dif- 
ferent filling fractions and studying excitation spectra in 
these states. However, in some of the systems of current 
theoretical and experimental interest, exact diagonaliza- 
tion has almost reached its limits because of the expo- 
nential growth of Hilbert-space dimension with the sys- 
tem size. There are many examples where a new method 
would be beneficial, such as QHE systems with spin de- 
grees of freedom and double-layer configurations, studies 
of higher filling fractions when several Landau levels are 
occupied as well as problems related to nonequilibrium 



dynamics of electrons in quantum Hall edge states which 
are still far from being well understood, see for example 
Ref.33. 

Recent advances in experiments with cold atomic gases 
allow to study systems of neutral particles, both fermions 
and bosons, under rotation, which is equivalent to hav- 
ing a coupling to a magnetic field in the rotating frame. 
These experiments showed formation of the triangu- 
lar Abrikosov vortex lattice in rotating Bose-Einstein 
condensates and in systems of fermions with attrac- 
tive interactions A fast rotation regime was achieved 
experimentallyji^ where disappearance of the vortex lat- 
tice was observed at very high rotation frequency. Al- 
ternative way of creating artificial magnetic fields based 
on the use of optical lattices was proposed in Rcf.^^ and 
has been recently realized in experimenlji^. An interest- 
ing question about melting of the vortex lattice under 
increased rotation was investigated by Cooper et ali^. 
Quantum Hall fractions of rotating bosons have been 
extensively studied with exact diagonalization by Reg- 
nault et ali^ Critical properties and the phase diagram 
of bosons in extremely elongated traps (narrow channel 
geometry) were discussed recently by Matveenko et alt^ 

In this paper we report on the DMRG technique for 
rotating bosons in the quantum Hall effect regime. This 
has been possible because of the special structure of 
the Hamiltonian projected onto the lowest Landau level 
which becomes essentially one-dimensional in the mo- 
mentum space with long-range interactions generated by 
the projection. It is well-known that DMRG works better 
in systems with short-range interactions, however there 
are successful applications of the method to calculations 
for strips of finite width. ^® Moreover in the case of quan- 
tum Hall effect in a cylinder geometry the matrix ele- 
ments of interactions between different momentum eigen- 
states are exponentially suppressed at large transferred 
momenta. We discuss convergence of the method in 
spherical geometry and on a cylinder for different filling 
fractions. 

The structure of the paper is the following. In Sect. 
II we introduce the Hamiltonian describing a system of 
fast-rotating bosons. In Sect. HI we explain our method 
on an example of the narrow-channel geometry. In Sect. 
IV we present the results of our calculations for systems 
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in spherical and cylinder geometries. 

II. BOSONS UNDER EXTREME ROTATION 

Let us consider a system of neutral bosonic atoms, 
which are strongly confined in the z direction by an 
external trapping potential with the frequency such 
that the bosons are in the ground state of the harmonic 
well. In this limit the system becomes essentially two- 
dimensional in the (x, y) plane. The bosons are confined 
in this plane by an additional harmonic trapping poten- 
tial F(r), with r = {a;, J/} which rotates around the z 
axis with a frequency = f2z. A single-particle Hamil- 
tonian for this systems is equivalent to one describing a 
particle with an effective charge q* and a mass m in an 
uniform magnetic field pointing along the z axis. In the 
case of a symmetric harmonic trap with a frequency uj 
and y(r) = muj'^Y^ /2 this Hamiltonian reads 

H^^{v-q*A/cf + \m{u:'-Q')v\ (1) 

Here p is a two-dimensional momentum operator in 
the {x,y) plane, A = {mc/q*)[^ x r] is a vector po- 
tential which corresponds to an effective magnetic field 
B = [V X A] = 2niVlc/q*7. and c is a speed of light. At 
the critical rotation frequency = w, the residual con- 
finement vanishes and the dynamics is governed by the 
Hamiltonian of a charged spinless particle moving in an 
infinite two-dimensional plane in external magnetic field. 

H,^^{v-q*A/cf. (2) 

In the following we will study systems of bosonic atoms 
with short-range delta-function interactions 

Htt\v-v')=g,a5^^''\v-v'\ (3) 

with g^d = iTrh^as/m, where Us is the three-dimensional 
scattering length. If t he harm onic oscillator length in 
the z direction, — y^h/muJz, is much larger than the 
scattering length \as\ and the characteristic radius of in- 
terparticle interaction, one can find for the effective in- 
teraction constanljiS 

The many-particle Hamiltonian of the system in the 
quasi-two-dimensional geometry reads 

H2d = ,^{Pa-q*Aa/c)' + Y,g2dS^"'\ra-n) (5) 
2m ^-^ 

a a<b 

When the chemical potential of the gas /i ~ g2dnp with 
the density rip, is much smaller than the cyclotron gap 
hujc = hq*B/mc = 2MI, the system can be effectively 
described by projecting it to the lowest Landau level. 



In this case the kinetic energy is quenched to zero and 
the dynamics is governed by interactions which leads to 
ground states depending on the boson filling fraction 
which is given by the ratio of the number of particles Np 
to the number of fiux quanta Ny — A/27rP, where A is 
the area of the system and I — ^Jhc/q*B = h/2mil is 
a magnetic length. In the limit of large filling fractions 
the ground state is gapless and is given by an Abrikosov 
vortex lattice. At small filling fractions it is represented 
by strongly correlated states which are gapped. The 
Abrikosov lattice and the gapped states are connected 
at intermediate filling fractions of the order of ~ 7 by 
a phase transition with melting of the vortex latticej^ 



III. DISCUSSION OF THE DMRG METHOD 

Density matrix renormalization group is a stan- 
dard method, which was initially developed for one- 
dimensional quantum lattice models by White in Refi^ 
It allows to calculate various quantum mechanical ob- 
servables such as energies, correlation functions, etc. for 
the ground and a few excited states of the systems with 
short-range interactions with an extremely high preci- 
sion. In the following we will extend the standard DMRG 
approach to study systems of rotating bosons described 
by the projected to the lowest Landau level Hamiltonian 
([S]). DMRG for fermions in the QHE regime have been 
studied in Refs.-*2i 

The general idea of the method is to use the eigenvec- 
tors, corresponding to the highest eigenvalues of the den- 
sity matrix, which is calculated for a part of the system, 
as an optimal basis to represent a target state, usually the 
ground state. There are two main versions of the algo- 
rithm i.e. the finite size and the infinite size methodsi^ 
To achieve best results one usually starts with the infinite 
size approach and correct it afterwards by implementing 
finite size sweeps. In our case because of the bosonic 
statistics particle occupation numbers can be large, more- 
over one has to impose constraints on the total number 
of particles and the total momentum, which makes it dif- 
ficult to use standard methods and we have to implement 
a different approach which is presented below. 



A. Narrow-channel geometry 

We will discuss our method using an example of the 
narrow-channel geometry, which was introduced^^ and 
studied recently-^ in the context of rotating bosons in 
strongly elongated traps in the mean-field regimeji^. This 
situation arises when the confining potential in the x di- 
rection is smaller than the trapping frequency in the y- 
direction and the rotation frequency is equal to uj^ ■ In 
this limit the system becomes infinitely elongated in the 
X direction and have remaining confinement in the y di- 
rection. The single-particle Hamiltonian in the rotating 
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frame in the Landau gauge read; 
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^0 = 7^(P + 2mnye.^f + l-muj'ty'^ , (6) 
2m 2 



■ w, 



2 - > 0. 



where up_ 

We impose periodic boundary conditions on the sys- 
tem of finite length L in the x direction so that the cor- 
responding momentum is quantized as fc = 21:71/ L with 
n G Z and assume that 57 ^ The eigenstates of 

Hamihonian (jS]) in the lowest Landau level are given by 
plane-waves with momentum k in the x direction and 
have a gaussian profile in the y direction which is cen- 
tered at the positions shifted by the value of kP with 
respect to the origin. These eigenfunctions read 
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(7) 



The set of states ([7]) with all possible momenta k repre- 
sents a basis on the lowest Landau level. 

The many-particle Hamiltonian for bosons with delta- 
function interactions ([3]) in the narrow channel geometry 
is given by the equation 



dxdy ^+ 



(p + 2mVlyex) 
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dxdy (8) 



In order to project this Hamiltonian onto the lowest Lan- 
dau level we write the operators ^{x, y) in the basis of 
eigenfunctions ^ 



(9) 



where ak,a'^ are boson annihilation and creation opera- 
tors with the commutation relations [a/j, a^] = Skq. After 
substitution of Eq. © into Eq. ^ and integration over 
x,y we obtain 

00 ^ 
Hnc= ^ k'^aldk + -'^Vijkatd1jakaik+],k+i- (10) 



ijkl 



Here the energy is given in units of fi^ /2m*P and — 
m{2Vl/uj^Y is the effective mass. The matrix elements 
of the interaction potential are given by the equation 



Vijk = ge 



-'-^[{^'kf + a-kf] 



(11) 



where g = 2 / TTlm* g2d / fi?' L is the dimensionless cou- 
pling constant. 

The properties of the ground states of the Hamilto- 
nian pop depend on three parameters for the system of 
finite size, namely the coupling constant the length 
of the system L and the total number of particles Np. 



a) IniliLil Wilson RG wamiup procudui 



b) Finile-sizc DMRG sweeps with one cenlral site 



FIG. 1: a) Warmup process using the standard Wilson renor- 
malization procedure. A sweep is started from the leftmost 
site with one site added at each step on the right until we 
reach the central site, b) Finite-size DMRG sweeps with a 
single central site after the initial Wilson RG procedure. 



In the thermodynamic limit the length dependence dis- 
appears and we will only have two parameters, the in- 
teraction strength and the dimensionless linear particle 
density Upl = Npl/L. The limit Upl <^ 1 corresponds 
to a quasi-one-dimensional case and is similar to a Lieb- 
Liniger gas, the opposite limit npl ^ 1 for small inter- 
actions (mean-field regime) was studied in i^^i^^ Here we 
will consider another regime of strong interactions g S> 1 
and large particle densities Upl ^ 1. 



DMRG algorithm for bosonic quantum Hall 
effect 



Our method is based on application of DMRG in the 
momentum space which was proposed by Xiang in the 
context of the Fermi-Hubbard modelSi. However, in our 
case the boson occupation numbers can take any value 
from zero to the total number of particles in contrast to 
fermions and one has to modify the method of Ref.™ In 
this subsection we will consider the system with Hamilto- 
nian (ITOl) with total number of particles Np and length L. 

We put the system on a lattice with a fixed num- 
ber of sites Ns = 2Nmax + 1 hi momentum space with 
k = 27:71/ where 71 takes values 71 = —Nmax--Nmax and 
Nmax is a cutoff. In general, the number of sites is infi- 
nite, but for the narrow-channel geometry there is always 
a physical cut-off due to the "kinetic energy" term in the 
Hamiltonian. The Hilbert space of the problem is rep- 
resented by the states with the sites filled with bosons 
taking occupation numbers 7ik from to Np. Hamilto- 
nian pU)) conserves total number of particles Np — J2 "-fc 
and total momentum -P = ^^fc j which we implement 
as constraints on allowed basis states. 

Step (I). We start with the Wilson renormalization 
group procedure from the leftmost site which has mo- 
mentum —Nmax- In the following we will use site num- 
bers to denote the momenta. We take few sites m with 
momenta from —Nmax to —Nmax + m — 1. In our cal- 
culations we used to from one to three. We diagonalize 
the Hamiltonian in the basis generated on these sites in 
all possible subspaces {71, p) of number of particles and 
momenta, where 71 is an integer which runs from zero 
to a total number of particles and p takes values from 
—n X Nmax to n X {—Nmax + 712 — 1). We reduce the 
basis by omitting the states which can not contribute to 
the target state. For example if we are interested in the 
state with total momentum P, we should not take into 
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account subspaces which have no partner {n',p') on the 
sites which do not belong to a current block, i.e. if there 
is no set of {n' ,p') such that n + n' = Np and p + p' — P- 
After diagonalizing the Hamiltonian in all the subspaces, 
we project it in every subspace together with the matrices 
corresponding to the matrix elements of the operators Op 
between allowed bra and ket subspaces, onto the lowest 
energy eigenstates of the Hamiltonian. In the calcula- 
tions we take up to four states in every subspace. After 
projection we save the operator matrices on a computer 
hard disk. 

Step (H). We add a single site to the block which is 
calculated in the previous step. When the first site is 
added it would be an initial block diagonalized in the Step 
I. This is done by generating new basis blocks {ni,pi) 
from the states of the previous block and the states on the 
site with particle number and momentum conservation 
constraints as in the Step I, here i is the current number 
of the Wilson RG step. We generate the Hamiltonian 
in the new combined block in all {ni,pi) subspaces by 
taking a tensor product of the operators, represented by 
matrices, from the previous block and from the added 
site. We project the operators and the Hamiltonian to 
its lowest energy eigenstates as in Step I and save them 
to the hard disk. In addition, every RG step we save the 
following set of operators 

Ao(i) = fli, j) = a+flj, (12) 

i2(*,j) = E«^+fc«»+'»^"''^'' (13) 

k 

A,{€) = Ea+afca.+,_,e-V2[(-'=)^+0-'=)^ (M) 

This procedure suggested by Xiang in Refj^ significantly 
reduces computation time. We repeat Step II until we 
reach the central site. It is important to mention that 
we keep all the matrix elements of the interactions in the 
Hamiltonian and allow for all possible boson occupation 
numbers, at a given total number of particles, in the basis 
states. 

Step (HI). In this step we start finite-size version of 
DMRG. First, we separate the Hamiltonian into two 
blocks and a single central site in the following way 

H = Hl+Hl,+H,+ H.r + Hr + Hlb., (15) 

where L, R denote the left and the right block, the site 
is represented by the symbol •. The site Hamiltonian is 
diagonal and is given by the sum of the kinetic energy 
term and the part of the interactions with operators act- 
ing on the site occupation numbers. The operators H^ r 
represent the contribution to the Hamiltonian from the 
left or the right block, the operators Hlr are composed 
from the contributions connecting the left and the right 
blocks, Hl* and HmR connect the central site to the left 
or the right block correspondingly. 

At every DMRG step we generate a new basis to ob- 
tain the set of states with the target quantum numbers 



{Np , P) from the possible states in the blocks which sat- 
isfy the following equations Np = N^ + N, + Nr and 
P = Pl+ P,+ Pr, where Nl^r,,, Pl.r,* are the number 
of particles and momentum in a given subspace of the 
corresponding block or site. 

The Hamiltonian for the right block is obtained by re- 
flecting the Hamiltonian of the left block. We create the 
full Hamiltonian of the system by a tensor product of the 
matrices from the left block, the site and the right block 
in the basis calculated above. We use the eigenstates 
of this Hamiltonian to calculate the density matrix for 
the left or right part of the system by tracing out the 
coordinates of the remaining part in every subspace of 
quantum numbers. After diagonalizing the density ma- 
trix we generate a matrix of the eigenvectors with highest 
eigenvalues adding the extra eigenvectors if necessary to 
have at least one eigenvector in every subspace. We fix 
the maximum number of eigenvectors or select the ones 
which eigenvalues are greater than some small number s. 
In our calculations we choose e in the range from 10~* to 
10^^ while checking convergence. The maximum dimen- 
sion of the Hamiltonian matrix which we diagonalized 
during DMRG sweeps was of the order of 10^. 

We also implemented single site corrections in the 
DMRG algorithm, which had been proposed by White in 
Refi^, by admixing information from the Hamiltonian to 
the density matrix. This is necessary because the inter- 
actions in the projected Hamiltonian are not short-range. 
These corrections improve the convergence by a consid- 
erable amount. Although in our calculations we used a 
version of the algorithm with a single central site, it is 
possible to extend it to two or more sites. However this 
would require more computer resources. 



IV. APPLICATIONS 

In this section we apply the method developed above 
to study systems of bosons with contact interactions on 
the lowest Landau level in cases of the narrow channel 
geometry in the limit of large densities and strong inter- 
actions as well as cylinder and a sphere. We discuss the 
convergence of the method for different filling fractions 
in various geometries. 



A. Bosons in the narrow channel geometry in the 
strongly interacting regime 

The system of bosons in the narrow channel is de- 
scribed by Hamiltonian (flUl) and was studied recently 
in Refs J^i^ in the mean-field regime in the limit of large 
density and weak interactions. It shows a series of phase 
transitions as a function of the interaction strength or the 
particle density, between the states with different num- 
ber of vortex rows. These states are gapless and have 
a Bogoliubov spectrum at small energies. In the oppo- 
site limit of very strong interactions bosons enter a state 
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FIG. 2: Density profiles for rotating bosons in tlie narrow 
channel geometry (in presence of a trap) for Np — 20, L = 10/, 
and different interaction strengths; long dashed line g — 100, 
dotted line g = 65 , short-dashed line g — 55, dot-dashed line 
<7 = 40. 



which is similar to the Laughlin v = 1/2 state, in the case 
of contact interaction potential, which minimizes their in- 
teraction energy. This state has a bulk gap and gapless 
edge states. When the interaction strength is decreased, 
the gap should close and there has to be be a phase transi- 
tion at intermediate values of interactions. Notice that in 
the narrow channel the system area is increased with in- 
creasing interaction strength as bosons repel each other, 
which leads to the increase in the number of flux quanta 
penetrating the system and decrease in the effective fill- 
ing fraction i/^ 

Here we calculate a few phases starting from very large 
interactions. The results for the integrated over x particle 
density n{y) , which is given by the equation 



n{y) 



1 



fe 



(16) 



where — (a^Ofc), are presented in Figl5] We take the 
system size L = 101, number of particles Np = 20 and 
up to Ns = 41 momentum states in our calculations. For 
very strong interactions g > 68 the ground state has a 
density profile which is very close to the one calculated for 
a cylinder geometry with v = 1/2 and which is described 
by the Laughlin wave-function 



(17) 



see FiglHl In this state the average occupation numbers 
rife ~ 1/2 in the bulk and increase at the edges by a small 
amount. The total number of occupied k-states is equal 
to Ns = 2Np — 1. When the interaction is decreased, 
the system enters another state through first order phase 
transition at ~ 68. In this state the system shrinks by 
exactly two momentum states ±kmax-, where k^ax is the 



state with the highest momentum which is occupied in 
the quasi-Laughlin wave-function, and the extra particles 
move to the center of the trap which is seen as a bump in 
the density on the Fig. [5]at g = 65. The area surrounding 
this bump continues to have the same constant density as 
before. This transition can be supported by the following 
argument. The energy cost of putting a particle to the 
center of the trap is of the order of g, and the gain of 
removing a particle from the state with momentum k^ax 
is of the order of A:^^,^. The transition happens when g ~ 
fc^jjj,. Decreasing the interactions further we observe a 
series of phase transitions with the analogous behavior 
This situation is similar to phase transitions in systems 
of lattice bosons in a trap, which shows wedding-cake 
structures which appear when changing the interaction 
strength or particle density. Similar structures have been 
studied by Cooper et. al.^^ in the case of quantum Hall 
effect. 



B. Cylinder geometry 

Let us now consider a system of bosons on a cylinder. 
This geometry was discussed in the fermionic case by 
Rezayi, Haldane'^^ and later studied by Bergholtz and 
Karlhede using DMRG in Ref.^^ 

We will impose periodic boundary conditions on the 
system of size L along the x axis, so that the momen- 
tum is quantized as fc = 27to/L where n G Z. The 
lowest Landau level basis wave-functions are given by 
the equation ([7]). This geometry could possibly be re- 
alized in experiments using proposals of creating artifi- 
cial magnetic fields in optical latticesJ^ The number of 
momentum eigenstates is finite because of the finite ex- 
tent W of the bar in the y direction and is given by 
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FIG. 3: Density profiles for the cylinder geometry with L = 
lO; for different filling fractions; i/ = 1/2 with Np = 20, Ns = 
39 solid line, i/ = 2/3 with Np — 22, A'^s ~ 31 dotted line and 
i/ = 1 with Np = 24, Na = 23 dot-dashed line. 
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FIG. 4: Density profiles for tlie cylinder geometry, A^s = 39, 
L = IQl and A'j, = 20 (dashed line) which corresponds to 
V = 1/2 Laughlin state, the same with one extra particle 
A'p = 21 (dotted line) and with two extra particles Np — 22 
(dot-dashed line). 



FIG. 6: Density profiles for the cylinder geometry with L = 
101 and Ns = 23 for A'p = 24 (dashed line) corresponding to 
1^ = 1 filling fraction, the same state with one particle removed 
A'p = 23 (dotted line) and one particle added A'p — 25 (dot- 
dashed line). 
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FIG. 5: Density profiles for the cylinder geometry with L — 
lOZ and Ns = 31 for A'p — 22 (dashed line) corresponding 
to = 2/3 filling fraction, the same state with one particle 
removed A'p = 21 (dotted line) and one particle added (dot- 
dashed line). 



system reads 



— ^ Vijk&fa'j'akaiSi+j^k+i, 



(18) 
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where matrix elements Vijk are given by the equation 
PT|) with g = g2d/V^Ll which plays a role of the energy 
scale. 

It is equivalent to the Hamiltonian of the narrow chan- 
nel geometry with "kinetic energy" term set to zero. 

We calculated the density profiles for systems of 
length L = 101 and different filling fractions, which is 
shown in Fig. [31 The energy of the i/ = 1/2 Laughlin 
state is zero, the best energies per particle of the v — 
2/3 state on a cylinder E^^2/3 ~ 0.489853(7 and for the 



1/ = 1 state -Ei/=i w 1.199865(7. These energies have 
been calculated with relative precision w 10~^, however 
convergence of the problem is nonlinear and the error is 
difficult to estimate. We have also calculated the energies 
and the density profiles for the states with extra added 
or removed particles which is presented in Figs. I4I5I6I 
For the Laughlin ly = 1/2 state the energy of the state 

with one extra particle is equal to Ae'^'^^^^ = 2.1917252g 
which corresponds to Np = 21 and for two particles is 
given by Asl't^-^'^^ = 3.5763744g for N, = 39 and L = 
101. For the = 2/3 state the energy cost of adding one 
par tide is A£;^+2/3 = 2.532379.g and the energy change 

for removing a particle Ai?^^2/3 — —1.5541128(7, here 
L = 101 and = 31. The energy change for the = 
1 state with added and removed particle are given by 



AE, 



(+1) 



^t.=i - 3.663195 and AeIJ' = -3.290275 with Np = 
25 and A'p = 23 particles correspondingly for Ng = 23 
and L = 101. 

In general the properties of the system will depend on 
the length L and to approach a thermodynamic limit one 
have to send L to infinity. In the case of = 1/2 finite- 
size effects in the density profile are already very small 
for L = 10, although one can still see oscillations for the 
states with higher filling fractions v = 2/3 and 1^=1. 
This could be an effect of the finite system length L or 
finite number of particles. In order to separate these two 
contributions it would be necessary to perform calcula- 
tions for systems of larger size. It is important to men- 
tion that with increasing L the interactions between par- 
ticles in different momentum states will increase which 
will result in decrease of the numerical convergence and 
will require more computer resources. However, it can 
be shown in certain cases that important physics can 
be inferred from studies of systems with finite lengths. 
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The authors of paper Ref.'^'' observed that a system of 
fermions on a cyUnder undergows continuous transforma- 
tion, i.e. without a phase transition, from the Laughhn 
to the Tao-Thoules state (see also Ref.^''). Similar results 
were obtained in Ref.^^ for v — Xjl state of electrons on 
a cyhnder and in more general case on a thin torus in 
Refj^ One of the interesting apphcations of our method 
would be to study this adiabaticity in the bosonic case. 



C. Spherical geometry 

Let us now turn to the case of the spherical geometry. 
Exact diagonalization on a sphere, which was introduced 
by Haldane in Ref is one of the most used in the finite- 
size studies of fractional quantum Hall effect for both 
bosons and fermions. Application of DMRG to systems 
of fermions on a sphere was developed in Refj^, here we 
consider bosonic case. 

The Hamiltonian of the system of bosons with delta- 
function interaction in the spherical geometry is given by 
the equation 



m4 f^mi +^TT'2 :"T'3+''^4 1 



(19) 



where the matrix elements of interactions on the lowest 
Landau level read 



, g (25 + i)MntiC|i^J^/^ ^ ^20) 



47rS' 45 + 1 



^4S 
^2S+mi+m2 



Here = ?i!/fc!(rt — h)\ are binomial coefficients and 
S = Nmax which gives A^s = 25 + 1. 



-1 



< 



^ -3 



Np=20, N=39, L=10, cylinder 
Np=16, Ns=31, sphere 
Np=16, N5=15, sphere 
Np=20, Ns=39, sphere 




50 

iteration # 



100 



FIG. 7: Absolute error in the ground state energy for the 
Laughlin state = 1/2 on a cylinder with L = lOi, Ap = 20, 
Aa = 39 (circles), v = 1/2 Laughlin state on a sphere for 
iVp = 16, = 31 (squares), Ap = 20, A^ = 39 (diamonds), 
1/ = 1 on a sphere for A'p = 16, As = 15 (triangles) as a 
function of the number of iterations. 



In this paper we studied convergence of the ground 
state energies at filling fractions v = \ and = 1/2 on a 
sphere. The results are presented in Fig. [7] together with 
the convergence for the v — Xjl Laughlin state in the 
cylinder geometry for comparison. In the calculations we 
used the states of the density matrix with the eigenval- 
ues larger than e — 10~*. We find that for a cylinder 
with h = lOZ, convergence is generally much better than 
that of the spherical geometry. This is because the ma- 
trix elements of the Hamiltonian on a cylinder for this 
value of L fall-off much faster with the distance between 
sites in the momentum space than for a sphere (see also 
discussion in Sec lIVBI) . 



CONCLUSIONS 



In conclusion, we developed a momentum space den- 
sity matrix renormalization group (DMRG) technique to 
study bosonic fractional quantum Hall effect. Using this 
method in the narrow channel geometry we found quan- 
tum phase transitions out of the = 1/2 Laughlin- type 
state at the critical value of the interaction strength. 
We calculated energies of ground states as well as states 
with added or removed particles and the density profiles 
at filling fractions v = 1/2, v = 2/3 and = 1 on a 
cylinder. The number of particles in our calculations is 
considerably larger than in current exact diagonalization 
studies of bosonic QHE. We found that convergence of 
the method strongly depend on geometry and it requires 
more computational efforts to study systems on a sphere, 
compared to a cylinder, provided that cylinder is thin 
enough (which is always the case in our calculations). 
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